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Abstract: The distribution of the magnitudes of damage avalanches during a failure process typically 
follows a power law. When these avalanches are recorded close to the point at which the 
system fails catastrophically, we find that the power law has an exponent which differs from 
the one characterizing the size distribution of all avalanches. We demonstrate this analytically 
for bundles of many fibers with statistically distributed breakdown thresholds for the individ- 
ual fibers. In this case the magnitude distribution D(A) for the avalanche size A follows a 
power law with ^ = 3/2 near complete failure, and £ = 5/2 elsewhere. We also study 
a network of electric fuses, and find numerically an exponent 2.0 near breakdown, and 3.0 
elsewhere. We propose that this crossover in the size distribution may be used as a signal for 
imminent system failure. 
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1 INTRODUCTION 

Catastrophic failures |H]|2j[3j|4j|5l are abundant in nature: earthquakes, landslides, mine-collapses, 
snow-avalanches etc. are well-known examples. A sudden catastrophic failure is a curse to hu- 
man society due to the devastation it causes in terms of properties and human lives. Therefore a 
fundamental challenge is to detect reliable precursors of such catastrophic events. This is also an 
important issue in strength considerations of materials as well as in construction engineering. 

During the failure process of composite materials under external stress, avalanches of differ- 
ent size are produced where an avalanche consists of simultaneous rupture of several elements. 
Such avalanches closely correspond to the bursts of acoustic emissions 10 [71 which are observed 
experimentally during the failure process of several materials. If one counts all avalanches till the 
complete failure, their size distribution is typically a power law. However we observed recently 
that if one records avalanches not for the entire failure process, but within a finite interval, a clear 
crossover behavior ||8] [9j is seen between two power laws, with a large exponent when the sys- 
tem is far away from the failure point and a much smaller exponent for avalanches in the vicinity 
of complete failure. Therefore such crossover behavior in the burst statistics can be taken as a 
criterion for imminent failure. In this report we discuss the crossover behavior in two different 
fracture-breakdown models- fiber bundle models [10-18] which describe the failure of composite 
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material under external load and fuse models |[T][T9j|20l which demonstrate the breakdown of elec- 
trical networks. An analytic derivation of the crossover behavior is given for fiber bundle models 
and numerical study confirms similar behavior in a fuse model. Thus we claim that such crossover 
behavior can be used as a signal of imminent failure. A recent observation [21 ] of the existence of 
a crossover behavior in the magnitude distribution of earthquakes within Japan has strengthened 
this claim. 



2 FIBER BUNDLE MODEL 



Our fiber bundle model consists of N elastic and parallel fibers, clamped at both ends, with statis- 
tically distributed thresholds for breakdown of individual fibers (Figure 1). The individual thresh- 
olds Xi are assumed to be independent random variables drawn from the same cumulative distri- 
bution function P(x) and a corresponding density function p(x): 



Prob(xj < x) = P(x) 



p(u) du. 



(1) 



Figure 1. A fiber bundle of N parallel fibers clamped at both ends. The externally applied force F corre- 
sponds to a stretching by an amount S. 



Whenever a fiber experiences a force equal to or greater than strength threshold X{, it breaks 
immediately and does not contribute to the strength of the bundle thereafter. The maximal load the 
bundle can resist before complete breakdown is called the critical load and its value depends upon 
the probability distribution of the thresholds. Two popular examples of threshold distributions are 
the uniform distribution 

_ , , | x lx r for < x < x r 

P{ x ) = S ■, r 

1 for x > x r , 

and the Weibull distribution 



(2) 



P(x) = 1 - exp(-(x/x r ) K ). (3) 
Here x r is a reference threshold, and the dimensionless number k is the Weibull index (Figure 2.) 
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Figure 2. The uniform distribution (A) and the Weibull distribution (B) with k = 5 (solid line) and k = 10 
(dotted line) . 

Fiber bundle models differ in the mechanism for how the extra stress caused by a fiber failure is 
redistributed among the unbroken fibers. The simplest models are the equal-load-sharing models, 
in which the load previously carried by a failed fiber is shared equally by all the remaining intact 
fibers in the system. In the present article we study the statistics of burst avalanches for equal- 
load-sharing models. 



2.1 Burst statistics 



The burst distribution D(A) is defined as the expected number of bursts in which A fibers break 
simultaneously when the bundle is stretched steadily until complete breakdown. Hemmer and 
Hansen performed a detail statistical analysis lfl4ll to find the burst distribution for this quasi-static 
situation. For a bundle of many fibers they calculated the probability of a burst of size A starting 
at fiber k with threshold value Xk as 



A 



A-l r 



A! 



1 



XkP(Xk) 

Q(xk) 



XkP(Xk) 



A-l 



x exp 



XkP{Xk) 



Q{xk) 



(4) 



Q{xk) 

where Q(x) = 1 — P(x). Since a burst of size A can occur at any point before complete break- 
down, the above expression has to be integrated over all possible values of xj., i.e, from to x c 
where x c is the maximum amount of stretching beyond which the bundle collapses completely. 
Therefore the burst distribution is given by lfT4ll 



D(A) _ A A -! e - A _. ri , A _, 



N A! 

where 



p(x)r(x)[l — r{x)\ exp [A r(x)] dx, (5) 



r(x) = 1 " ~Q{x) = Q{x)d-x [xQ{x)] - (6) 
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The integration yields the asymptotic behavior 



D(A)/N oc A~2, (7) 

which is universal under mild restrictions on the threshold distributions lfl5ll . As a check, we have 
done simulation experiments for different threshold distributions. Figure 3 show results for the 
uniform threshold distribution (0 and the Weibull distribution ® with index k = 5. 
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Figure 3. The burst distribution D(A)/N for the uniform distribution (A) and the Weibull distribution with 
index 5 (B). The dotted lines represent the power law with exponent —5/2. Both figures are based on 20000 
samples of bundles each with N = 10 6 fibers. 

2.2 Crossover behavior near failure point 




A 



Figure 4. The distribution of bursts for thresholds uniformly distributed in an interval (xo, x c ), with xq = 
and with xq = 0.9x c . The figure is based on 50 000 samples, each with N = 10 6 fibers. 

When all the bursts are recorded for the entire failure process, we have seen that the burst 
distribution D(A) follows the asymptotic power law D oc A~ 5 / 2 . If we just sample bursts that 
occur near criticality, a different behavior is seen. As an illustration we consider the uniform 
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threshold distribution, and compare the complete burst distribution with what one gets when one 
samples merely burst from breaking fibers in the threshold interval (0.9x c ,x c ). Figure 4 shows 
clearly that in the latter case a different power law is seen. 

From Eq. [6] we see that r(x) vanishes at the point x c . If we have a situation in which the 
weakest fiber has its threshold xq just a little below the critical value x c , the contribution to the 
integral in the expression © for the burst distribution will come from a small neighborhood of 
x c . Since r(x) vanishes at x c , it is small here, and we may in this narrow interval approximate the 
A-dependent factors in ([5]> as follows 



(1-rfe 



A Ar 



exp [A(ln(l — r) + r)] 

exp[-A(r 2 /2 + C(r 3 ))] » exp \-Ar(x) 2 /2 



We also have 



r(x) sa r'(x c )(x - x c ). 
Inserting everything into Eq. ((5]>, we obtain to dominating order 

D(A) 



(8) 
(9) 



N 



A A-1 e~ A 
A! 

A A-2 e -A p ( a . c ) 

\r'(x c )\Al 
A A - 2 e- A p{x c 



p(x c ) r'{x c )(x - x c )e- Ar ' {Xc) {x ~ Xc) /2 dx 



A! 



\r'(x c )\ 



-Ar'(x c ) 2 (x-x c ) 2 /2 



1 - e " A / A < 



with 



A, 



(10) 



(11) 



r'(x c ) 2 (x c - xq) 2 ' 

By use of the Stirling approximation A! ~ A A e~ A \/2-7rA, the burst distribution (fTOl) may be 

D(A) 



written as 



N 



CA- 5 ' 2 (l - e" A / Ac ) , 



with a nonzero constant 



C = (27T)- 1 / 2 p(x c )/\r'(x c )\. 
We can see from (fT2l) that there is a crossover at a burst length around A c : 

D(A) f A" 3 / 2 forA<A c 
TV ^ I A" 5 / 2 for A > A c 



(12) 
(13) 

(14) 



We have thus shown the existence of a crossover from the generic asymptotic behavior D oc 
A~ 5 ' 2 to the power law D oc A~ 3 / 2 near criticality, i.e., near global breakdown. The crossover is 
a universal phenomenon, independent of the threshold distribution p(x). 
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Figure 5. The distribution of bursts for the uniform threshold distribution (left) with xq = 0.80x c and for a 
Weibull distribution (right) with xq = 1 (square) and xq — 1.7 (circle). Both the figures are based on 50000 
samples with N — 10 6 fibers each. The straight lines represent two different power laws, and the arrows 
locate the crossover points A c ~ 12.5 and A c ~ 14.6, respectively. 

For the uniform distribution A c = (1 — xq/x c )~ 2 /2, so for xo = 0.8 x c , we have A c = 12.5. 
For the Weibull distribution P(x) = 1— exp(— (x— l) 10 ), where 1 < x < oo,wegetx c = 1.72858 
and for xq = 1.7, the crossover point will be at A c ~ 14.6. Such crossover is clearly observed 
(Figure 5) near the expected values A = A c = 12.5 and A = A c = 14.6, respectively, for the 
above distributions. 
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Figure 6. The distribution of bursts for the uniform threshold distribution for a single fiber bundle with 
N = 10 7 fibers. Results with xq = 0, i.e., when all avalanches are recorded, are shown as squares and data 
for avalanches near the critical point (xq = Q.9x c ) are shown by circles. 

The simulation results we have shown so far are based on averaging over a large number of 
samples. For applications it is important that crossover signal can be seen also in a single sample. 
We show in Figure 6 that equally clear crossover behavior is seen in a single fiber bundle when 
iV is large enough. Also, as a practical tool one must sample finite intervals (x«, Xf) during the 
fracture process. The crossover will be observed when the interval is close to the failure point (9|. 
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3 BURST STATISTICS IN THE FUSE MODEL 
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Figure 7. A fuse model and the burst distribution: System size is 100 x 100 and averages are taken for 300 samples. 
On the average, catastrophic failure sets in after 2097 fuses have blown. The circles denote the burst distribution 
measured throughout the entire breakdown process. The squares denote the burst distribution based on bursts 
appearing after the first 1000 fuses have blown. The triangles denote the burst distribution after 2090 fuses have 
blown. The two straight lines indicate power laws with exponents £ = 3 and 5 = 2, respectively. 

To test the crossover phenomenon in a more complex situation than for fiber bundles, we have 
studied burst distributions in the fuse model [1]. It consists of a lattice in which each bond is a 
fuse, i.e., an ohmic resistor as long as the electric current it carries is below a threshold value. If 
the threshold is exceeded, the fuse burns out irreversibly. The threshold t of each bond is drawn 
from an uncorrected distribution p(t). The lattice is placed between electrical bus bars and an 
increasing current is passed through it. The lattice is a two-dimensional square one placed at 45° 
with regards to the bus bars, and the Kirchhoff equations are solved numerically at each node 
assuming that all fuses have the same resistance. 

When one records all the bursts, the distribution follows a power law D(A) oc A~^ with £ = 3, 
which is consistent with the value reported in recent studies |fi~9ll20ll . We show the histogram in 
Figure 7. With a system size of 100 x 100, 2097 fuses blow on the average before catastrophic 
failure sets in. When measuring the burst distribution only after the first 2090 fuses have blown, a 
different power law is found, this time with £ = 2. After 1000 blown fuses, on the other hand, £ 
remains the same as for the histogram recording the entire failure process (Figure 7). 

In Figure 8, we show the power dissipation E in the network as a function of the number of 
blown fuses and as a function of the total current. The dissipation is given as the product of the 
voltage drop across the network V times the total current that flows through it. The breakdown 
process starts by following the lower curve, and follows the upper curve returning to the origin. 
It is interesting to note the linearity of the unstable branch of this curve. In Figure 9, we record 
the avalanche distribution for power dissipation, Dd(A). Recording, as before, the avalanche 
distribution throughout the entire process as well as recording only close to the point at which 
the system catastrophically fails, result in two power laws, with exponents £ = 2.7 and £ = 1.9, 
respectively. It is interesting to note that in this case there is not a difference of unity between the 
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two exponents. The power dissipation in the fuse model corresponds to the stored elastic energy 
in a network of elastic elements. Hence, the power dissipation avalanche histogram would in the 
mechanical system correspond to the released energy. Such a mechanical system could serve as a 
simple model for earthquakes. 
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Figure 8. Power dissipation £ as a function of the number of broken bonds (left) and as a function of the 
total current I flowing in the fuse model (right). 

The Gutenberg-Richter law (T] |3 relating the frequency of earthquakes with their magnitude 
is essentially a measure of the elastic energy released in the earth's crust, as the magnitude of an 
earthquake is the logarithm of the elastic energy released. Hence, the power dissipation avalanche 
histogram Dd(A) in the fuse model corresponds to the quantity that the Gutenberg-Richter law 
addresses in seismology. Furthermore, the power law character of D d(A) is consistent with the 
form of the Gutenberg-Richter law. It is then intriguing that there is a change in exponent £ also 
for this quantity when failure is imminent. 




Figure 9. The power dissipation avalanche histogram D^(A) for the fuse model. The slopes of the two 
straight lines are —2.7 and —1.9, respectively. The circles show the histogram of avalanches recorded 
through the entire process, whereas the squares show the histogram recorded after 2090 fuses have blown. 
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4 CONCLUDING REMARKS 



Establishing a signature of imminent failure is the principal objective in our study of different 
breakdown phenomena. The same goal is of course central to earthquake prediction scheme. As 
in different failure situations, bursts can be recorded from outside -without disturbing the ongoing 
failure process, burst statistics are much easily available and also contain reliable information of 
the failure process. Therefore, any signature in burst statistics that can warn of imminent system 
failure would be very useful in the sense of wide scope of applicability. The crossover behavior 
in burst distributions, we found in the fiber bundle models, is such a signature which signals that 
catastrophic failure is imminent. Similar crossover behavior is also seen in the burst distribution 
and energy distributions of the fuse model. Most important is that this crossover signal does not 
hinge on observing rare events and is seen also in a single system. Therefore, such signature has 
a strong potential to be used as useful detection tool. It should be mentioned that most recently, 
Kawamura ||2T1 has observed a change in exponent values of the local magnitude distributions of 
earthquakes in Japan, before the onset of a mainshock (Figure 10). This observation has definitely 
strengthened our claim of using crossover signals in burst statistics as a criterion for imminent 
failure. 
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Figure 10. Crossover signature in the local magnitude distributions of earthquakes in Japan . The exponent 
of the distribution during 100 days before a mainshock is about 0.60, much smaller than the average value 
0.88 ED- 
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